Erosion waves: transverse instabilities and fingering. 
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Two laboratory scale experiments of dry and under-water avalanches of non-cohesive granular 
materials are investigated. We trigger solitary waves and study the conditions under which the 
front is transversally stable. We show the existence of a linear instability followed by a coarsening 
dynamics and finally the onset of a fingering pattern. Due to the different operating conditions, both 
experiments strongly differ by the spatial and time scales involved. Nevertheless, the quantitative 
agreement between the stability diagram, the wavelengths selected and the avalanche morphology 
reveals a common scenario for an erosion/deposition process. 

PACS numbers: 47, 47.10.+g, 68.08.-p, 68.08.Bc 



Avalanching processes leading to catastrophic trans- 
port of various natural materials do not only occur in the 
air as we know of snow avalanches, mud flows and their 
catastrophic human and economical toll. Such events 
frequently happen below the see level as they take many 
forms from turbidity currents to thick sediment waves 
sliding down the continental shelf. This is a fundamen- 
tal feature shaping the submarine morphology. From the 
modeling of risks point of view, important questions still 
remain such as to evaluate to which extend an initial trig- 
gering event (an earth quake, an eruption..) would be re- 
sponsible for a subsequent process that might propagate 
or amplify over large distances as an unstable matter 
wave. Unfortunately, the dynamics of such catastrophic 
events remains an issue so far lacking of conceptual clar- 
ity Q Q since (i) the rheology of the flows involved in an 
avalanche is complex and still unraveled, (ii) the physics 
of erosion/deposition mechanisms is essentially limited to 
empirical descriptions based on dimensional analysis and 
semi-empirical formulations. 

Recently, extensive laboratory-scale experiments on 
dry granular materials have aimed to unify the rheol- 
ogy of dense particulate flows in different geometries 
J3j. These flows can be organized into two sub-classes, 
quasi-static creep flows characterised by an exponential- 
like velocity profile, and shear flows, whose rheology is 
local 0, [H El- In the case of fully developed granu- 
lar flows, measurements are now sufficiently precise and 
reproducible to evidence dependencies of the constitu- 
tive relation on microscopic granular features like rough 
sandy material vs round spherical particles 0, @ ■ But 
still, a full understanding of the flow constitutive rela- 
tions would require a pertinent description of the passage 
between the blocked and the mobile phases (the jamming 
transition). For instance, the onset and the dynamics of 
triangular shape avalanches was shown to be controlled 
essentially by the presence of a metastable substrate and 
the interplay between erosion and deposition processes 
0, ^3 • From a theoretical point of view, descriptions of 
the granular bed mobilisation have been proposed, either 
empirically ('a la Saint-Venant') |lllll2| or phenomeno- 
logically by phase field methods [l£ 
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FIG. 1: a Experimental set-up. b Stability diagram: h at0 p 
is the thickness of the sediment left after an avalanche for 
a given angle 8, in air (•) and in water (O); h s tart(0) is 
the maximum stable height of sediment, in air (T) and in 
water (v). h 3ta rt(0) and h s t op (0) are fitted by the form 
h — 61n((tan# — fi)/8fi) (solid lines). Avalanches triggered 
in region II are stable while they exhibit a transverse insta- 
bility in region III. In particular, solitary erosion waves are 
evidenced when starting from the stable height h 3 t op - c Front 
profile x(y) obtained after image processing by a correlation 
technique, d The corresponding correlation function C(y) al- 
lows to define the average wavelength A and amplitude A. 



Avalanche fronts flowing on solid rough substrates are 
transversally stable, the transverse coupling due to grav- 
ity being essentially a stabilizing mechanism [Til fl5| . 
But, in the particular case of a strong size bidispersity, 
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FIG. 2: Solitary erosion wave profile Sh = h—h stov rescaled by 
d (dotted line) and surface velocity profile v rescaled by \fgd 
(solid line) (dry, 6 — 32°, region II). Inset: spatio-temporal 
diagram done with a fast camera (125 Hz), showing the parti- 
cle motion as well as the profile height (deflection of the laser 
sheet). It can be observed that the surface grain velocity 
tends at the front towards the solitary wave velocity v a 

an avalanche front may exhibit a fingering pattern ex- 
plained by a pinning mechanism 16, 17] : the grains of 
larger size gathering in the finger thrusts are hindering 
locally the avalanche front progression. In this paper, we 
present an experimental study of avalanche fronts devel- 
oping over an erodible granular substrate, in the air and 
under water. We demonstrate the existence of a linear 
transverse instability of solitary erosion waves, although 
the rough grains we use exhibit a narrow polydispersity. 

The avalanching set-ups consist of a thin layer of grains 
deposited on a substrate that can be tilted at a value 9. 
The dry gran ular set-up is similar to the one of Daerr 
et al. [3, Ey]- The avalanche track is 70 cm wide and 
120 cm long. The granular medium is Fontainebleau sand 
of a medium size d = 300 /im and the track bottom is 
made of black velvet. For under water avalanches, the 
set up size is quite smaller. The avalanche track is the 
bottom of a plexiglass tank that can be tilted up to an 
angle 9 from an horizontal position. The avalanche track 
width is 15 cm and so is the track length. The granu- 
lar sediment is an aluminum oxide powder of size either 
d = 30 fim or 40 /im. To avoid interparticle cohesion, it 
is sufficient to maintain the pH value close to 4 by ad- 
equate addition of hydrochloric acid ^| . The substrate 
is initially set at an horizontal position and a fixed mass 
of powder is poured and suspended by vigorous stirring. 
A uniform sediment layer of height h then forms within 
10 min. The bottom is an abraded but transparent plex- 
iglass plate which offers the possibility to monitor the 
avalanche dynamics by transparency when illuminated 
from below. The profile of the avalanche front h(x,t) is 
obtained with a laser slicing technique and is resolved 




FIG. 3: Flowing part of solitary waves visualised by image 
difference (air, d = 300 /jm, 8 — 35 °, time interval 1.1 s), 
starting from a flat bed (left) or from an initial bed presenting 
a forced wavelength A = 6.5 cm. 

within 30 /im (0.1 d). The front dynamics is quantita- 
tively monitored by image processing of the avalanche 
front pictures. The front line equation x{y^) is then 
extracted (fig. Ut) an d the front line auto correlation 
function C(y,t) =< x(y + y',t)x(y' ,t) > y i is computed. 
Then, the correlation function first maximum is identified 
from which wc define the average wavelength A and the 
amplitude A = 2- s /2C(X) (fig. El). In addition, for dry 
avalanches, we measure the surface velocity field using a 
Particle Image Velocimetry technique. 

It has been shown that the stability of dry granular 
layers of depth h lying on a substrate inclined at an an- 
gle 9 can be simply apprehended by a diagram with two 
branches (fig-Gb) h s t a rt(9) and h s top(9) with the fol- 
lowing interpretation: a uniform deposit of height h will 
globally loose stability if tilted above the angle 9 defined 
by h = h sta rt{9) and the avalanching process will leave 
at rest a deposit of height h stop {9). The h start and h stop 
curves diverge at an asymptotic angle limit, respectively 
equal to the avalanche angle of the granular pile 9 a and 
to the repose angle 9 r . Between the two, a domain of 
metastability for the granular deposit is present. Inter- 
estingly, the stability curves obtained for dry and under- 
water layers bear the same features and fall on the same 
curve when the deposited height is rescaled by the grain 
size (fig. nj)). 

To initiate avalanche fronts both in air and under wa- 
ter, we designed a 'bulldozer' technique where a plate per- 
pendicular to the avalanche track scrapes the sediment 
at a constant velocity (fig. QJO . Although our results on 
avalanche stability are valid in the whole metastable re- 
gion (fig. E>) , we will limit ourselves here to experiments 
started from a stable sediment layer of height h stop (9). 
Once an autonomous avalanche front separates from the 
plate, the bulldozer driving stops. We obtain both for 
dry and under water avalanches the following results. 

(i) For 9 < 9 r the sediment is quite stable as an 
avalanche front cannot propagate autonomously down 
the slope (region / on graph lb): the perturbation is 
bound to fade away when the driving stops. 

(ii) For 9 r < 9 < 9 a , we always obtain transversely sta- 
ble autonomous avalanche fronts. We observed that the 
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FIG. 4: Time evolution of the wavelength A (water, d — 
40 pm, 9 = 37.1°) in a single typical realisation (•) and 
averaged over realisations (solid line) - the shadow zone in- 
dicates the standard deviation. After a small plateau at the 
initial wavelength Ao, A increases due to merging processes 
(lower photograph) until the value A^ which corresponds to 
the formation of non-interacting fingers (upper photograph). 




FIG. 5: Time evolution of the amplitude (air, d = 300 /im, 
9 = 35°) when the initial condition is forced at a given wave- 
length A = 12 mm (O), A = 30 mm (□), A = 90 mm (A) 
and A = 178 mm (T). Inset: linear growth rate a as a func- 
tion of the wave number k. The solid line is the best fit by 
a — cr m |fe|Ao(l — (fcAo) 2 /3), with a maximum growth rate 
a m ~ 2.5 s _1 for Ao — 4 cm. 



avalanche quickly converges toward a form which then re- 
mains constant. Furthermore, this solitary wave is found 
to be quite insensitive to the avalanche preparation de- 
tails within a range of scraping velocities or initial masses 
set into motion. For this systematic study, we have kept 
a constant scraping velocity at about one-third of the 
typical avalanche velocity v a . For each value of the - 
unique - control parameter 8, there is thus a single pos- 
sible solitary erosion wave. In water, v a is of the order 

• A/9 gd 2 , Ap 

of the Stokes velocity ~ 2 mm/s where = 3 

Pw 18 v p w 
is the density contrast between grains and water, v the 
water kinematic viscosity and g the gravity acceleration. 
In the air, the propagation velocity is of the order of 
\fgd ~ 5 cm/s. In figure [21 we show local sediment 
height h and local surface velocity v profiles for such an 
avalanche. It is important to note that the relation be- 
tween h and v does not correspond to the curve found for 
homogeneous steady flows 0- We rather obtain a rela- 
tion of the form v oc y/g(h — h s top)- This result suggests 
that the grains are not flowing down to the bottom and 
consequently, erosion/deposition processes are determi- 
nant to select the avalanche dynamics and shape. 

(Hi) For 8 > 6 a the neutral wave fronts are transver- 
sally unstable. After the initial instability, we have iden- 
tified a sequence of fusion processes increasing the spatial 
modulation lengths (coarsening scenario). Finally, the 
transverse destabilization ends up as a fingering pattern. 
In this final stage, the flowing zones are disconnected one 
from the others so that the wavelength does not evolve 
anymore. On figure^] we display a typical time evolution 
of the dominant wavelength extracted from the correla- 
tion function. In inset, a typical fusion event is displayed 
to illustrate the coarsening scenario. Because of the com- 
petition between unstable modes and the coarsening pro- 



cess taking place, the identification of a generic scenario 
for the transverse instability is problematic. This is the 
reason why, in addition to the experiments started from a 
flat bed we just described, we performed series of exper- 
iments starting from a modulated initial condition. The 
modulation at a given wavelength is simply produced by 
imprinting on the sediment surface regularly spaced thin 
scarifications. We find that the forced modes always fade 
away in region //, but on the other hand, in region 
the front modulations amplifies exponentially for a wide 
band of modes. The linear regime is clearly evidenced 
over one decade in amplitude. Non-linear effects start 
being visible when the amplitude becomes centimetric. 
The inset of figure [5] shows the dispersion relationship 
deduced from these experiments, which demonstrates the 
existence of an initial long wavelength linear instability. 

For experiments both in the air and under water per- 
formed in the unstable regime, we extract the two char- 
acteristic wavelengths. The initial wave length Ao would 
correspond, to the best of our experimental possibilities, 
to the fastest growing mode of the linear regime. Then, 
the wave length Aoo is taken at the onset of the finger- 
ing instability. In fig. we display both wavelengths 
rescaled by the grains sizes : Ao/d and Xoo/d, as a func- 
tion of the inclination angle 8. The selected wavelengths 
are typically larger than a grain size by at least two or- 
ders of magnitude. Note that the largest wavelengths 
measured are of the order of the track width (1800 d in 
water and 750 d in air). Furthermore, in the limit of fi- 
nite size effects and measurements uncertainties, we find 
that a value 8 = 8 a corresponds to a diverging bound- 
ary for the initial wavelength Xo/d. It suggests that Ao 
could directly scale on h start (solid line in fig. . Another 
striking feature is the collapse, on the same curve, of all 
data both in the air, underwater and for grains of differ- 
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FIG. 6: Initial (•) and final (O) wavelengths rescaled by d 
as a function of 6 (d — 40 /im, in water). In air, the wave- 
length obtained from the growth rate maximum (■) coincides 
with the initial wavelength measured in water (•). The final 
wavelength data in air (□) and water (O) also superpose. 
The error bars correspond to the dispersion of the data from 
a realisation to the other. As Ao diverges at 9 a , we have su- 
perimposed the curve 10 h sta rt(d) (solid line), which is a good 
approximation of Ao to the first order. The dotted line is the 
best fit of the final wavelength Aoo by the same logarithmic 
form as h sta rt(0) or h B t op (0). 

ent size. In the range of parameters where the fingering 
regime is reached before the end of the track, the ra- 
tio of the final to the initial wavelength is approximately 
constant and equal to Aoo/Ao ~ 3.5. The presence of a 
fingering instability is a quite fascinating feature of this 
avalanching process. Here, the fingering front stems from 
the onset of localized propagating waves following the 
transverse instability regime. These fingers are localized 
matter droplets with levees on the side and propagating 
in a quasi solitary mode and when they are fully devel- 



oped, their selected width is found to be quite sensitive 
to the slope (approximatively Ao, for both underwater 
and dry situations). 

In this letter, we have investigated the dynamics of 
underwater and dry granular avalanches taking place on 
a erodible substrate. We have identified the domain of 
existence for solitary waves going down the slope with- 
out changing form. For angles larger than the avalanche 
angle, we proved the existence of a linear transverse in- 
stability which further develops via a coarsening fusion 
process and finally ends up as a fingering pattern. The 
existence of solitary waves whose properties depend on 
the angle provides a new important test to models. For 
instance, it is easy to show that they cannot be recovered 
in Saint- Venant models like [l4| which do not include a 
static erodible layer below the avalanche. The relation- 
ship between h and v provides an indirect evidence that 
a static layer does exist close to the front and below the 
tail of the avalanche. The mechanism responsible for the 
instability remains yet to be identified. Nevertheless, the 
inhibition of this instability on a solid bottom as well as 
the scaling of the initial wavelength Ao on h s t a rt suggests 
that erosion/deposition processes in the avalanche depth 
play a determinant role. Solving the problem experi- 
mentally is difficult as it would require an experimental 
access to the jammed and rolling heights in order to de- 
termine transverse modulations. In the final stage of the 
instability, fingers appear as droplet like solutions of the 
erosion/deposition process and thus look essentially dif- 
ferent from the segregation fingers reported on a rough 
substrate ^(|. Note that their shape is reminiscent of 
many natural patterns obtained in debris or mud flows 
[l7| which also display surprinsigly well selected widths 
at values about hundreds of a typical rock size. 

We thank P. Claudin, O. Pouliquen, A. Daerr and 
I. Aranson for many discussions and suggestions on the 
experiment. 
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